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Abstract 

We study the general model of self-financing trading strategies in illiquid markets 
introduced by Schonbucher and Wilmott, 2000. A hedging strategy in the framework 
of this model satisfies a nonlinear partial differential equation (PDE) which contains 
some function g(cc). This function is deep connected to an utility function. 

We describe the Lie symmetry algebra of this PDE and provide a complete set of 
reductions of the PDE to ordinary differential equations (ODEs). In addition we are 
able to describe all types of functions g(a) for which the PDE admits an extended 
Lie group. Two of three special type functions lead to models introduced before by 
different authors, one is new. We clarify the connection between these three special 
models and the general model for trading strategies in illiquid markets. We study with 
the Lie group analysis the new special case of the PDE describing the self-financing 
strategies. In both, the general model and the new special model, we provide the 
optimal systems of subalgebras and study the complete set of reductions of the PDEs 
to different ODEs. In all cases we are able to provide explicit solutions to the new 
special model. In one of the cases the solutions describe power derivative products. 
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1 Introduction 

We study self-financing hedging strategies in the framework of a reaction-function 
model. The main subject of this model is a smooth reaction function \ff that gives the 
equilibrium stock price S t at time t as function of some fundamental value F, and the 
stock position of a large trader. In the framework of this model there are two types 
of traders in the market: ordinary investors and a large investor. The overall supply 
of the stock is normalized to one. The normalized stock demand of the ordinary 
investors at time t is modeled as a function D(F t ,S t ), where S t is the proposed price 
of the stock. The normalized stock demand of the large investor is written in the form 
P&t\ p > is a parameter that measures the size of the trader's position relative to the 
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total supply of the stock. The equilibrium price S t is then determined by the market 
clearing condition 

D(F„S t )+p* t = l. (l.D 

The equation ( |1.1[ ) admits a unique solution S t under suitable assumptions on the 
function D(F t ,S t ). Hence S t can be expressed as a function y of F t and pO,, so that 
5 f = Y(F„p& t ). 

Now we turn to the characterization of self-financing hedging strategies in the 
framework of the reaction-function model. Throughout we assume that the funda- 
mental value process F t follows a geometric Brownian motion with volatility a as in 
the Black-Scholes model. Moreover, we assume that the reaction function y{F t , p<& t ) 
is of the form 

V(f,a)=fg(a) (1.2) 
with some increasing function g(a). This holds for any model where 

D(f,s) = U(f/s) (1.3) 

for a strictly increasing function U : (0,°°) — > E with a suitable range. 

Assuming as before that the normalized trading strategy of the large trader is of 
the form p<&(t,S) for a smooth function <!>, we get from Ito's formula that 

dS t =g(p<f>(t,S t ))dF t + pF tga (p<l>(t,S t ))<Z> s (t,S t )dS t + b(t,S t )dt (1.4) 

(since S t = g(p<&(t,S t ))F t ). The precise form of b(t,S t ) is irrelevant to our purposes. 
Assume now that 

(l-pF t g a (p&(t,S t ))*sM))>0 a.s.. (1.5) 

This can be viewed as an upper bound on the permissible variations of the large 
trader's strategy. Rearrangement and integration (l — pF t g a (p<t>(t,St))<t>s(t,S t )) 1 
over both sides of equation ( 1.4 1 gives us the following dynamics of S: 

dSt = 1 -aS t dW t +b(t,S t )dt, (1.6) 

again the precise form of b is irrelevant. 

If we apply the Ito-Wentzell formula [ 1 ] to this equation we obtain the following 
PDE for a self-financing hedging strategy 



Under some special choices of the function g(oc) equation (1.7) represents the 
models introduced earlier. 



2 



If we takeg(a) = C2e C[(X , c\,c% — const, we obtain the model introduced in [ 10]. This 
model was studied with Lie group methods in Q, ||8l and some generalization of this 
model in J6J. 

In the case of g(oc) = C2(p + kcc) c i , c\ ,C2, k — const, equation coincides with 
the models introduced in [ 15 ] and studied with Lie group methods in [4] and Q. 

The main goal of this paper is to study the analytical properties of equation ( |1.7[ ) 
in the first line with the Lie group method. We will clarify which symmetry proper- 
ties has this equation and under which additional conditions on the function g(tt) it 
admits a richer symmetry group. 

If the symmetry group admitted by an equation is known then it can be used in 
many cases to reduce the given equation to a simpler one, for instance, to an ODE 
or even solve it. Solutions of such reduced equations are called invariant solutions 
because they are invariant under the action of the corresponding subgroup of the 
symmetry group of this equation. We follow this strategy and describe the complete 
set of reductions for the general model ( |1.7| ) in section [3] 

In addition we obtain also some special forms of utility functions for which the 
corresponding equation admits a richer symmetry group. The connection between the 
function g(cc) and the utility function is described in section [2j For a special form of 
utility function we obtain a special model for hedging strategies in illiquid markets. 
This special model is studied in section [4] Also for this model we present an optimal 
system of subalgebras, the complete set of reductions to ODEs. In this case we are 
able to solve all the ODEs and provide explicit formulas or graphs of solutions. 



2 The function g(ot), the corresponding utility function and 
the self-financing hedging strategy 

We use the classical assumptions (see ||9l) on a utility function U(x), modeling the 
utility of a market participant's wealth x at fixed time. 

Definition 2.1 Let U (x) is a utility function. Then U : R — > RU {— °°} 

1. is increasing on R, continuous on {U > —°°}, differentiable and strictly con- 
cave on the interior of {U > — °°}, 

2. U'{°°) =\im x ^ 00 U'{x) = 0, 

3. if a negative wealth is not allowed then U (x) = — °°,/orjc < 0, U (x) > —oo^for 
x > 0, and the Inada condition U'{0) = lim x ^o,x>o U'(x) = °° holds, 

4. if a negative wealth is allowed then U(x) > — °°, for i£l and U'(—°°) = 
]im x -+- OB U' (x) = °° holds. 
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Using the form of the factor representation of the reaction function (1.2) and ( 1.3 1 we 
obtain the following equation for the utility function 



U[— =l-a. (2.8) 

We restrict our consideration further to the case where a negative wealth is not 
allowed and consider the utility function U as a map U : M + ->lu{-t»} with the 
assumptions ([T]) - Q listed in the Definition 2.1 above. 



Because of this we can restrict the argument of the utility function U (x) on 



hence we look for a function g(tt) > 0. We are able to rewrite the relation (2.8) in 
a more convenient and explicit form for some typical special forms of the function 

g(a): 

1. If the function g(tt) takes the form g(tt) = C2e C[CX then it corresponds to the 
utility function which is equal to 

g(a)=c 2 e c,a , -> U(x) = -ln;t+ Cl+lnC2 , c h c 2 ^0. (2.9) 

Cl Cl 

Because the argument of the utility function should be positive the function 
g(a) = C2e c>a should be positive too, i.e. the coefficient C2 should be strictly 
positive C2 > 0. According to the assumption ([T]) of Definition 2.1 the utility 



u t + -- z*—2=0- (2.10) 



function U (x) should be an increasing function, then ci > 0. It is easy to prove 
that all other assumptions are as well satisfied for the utility function of the 
form $L9\ by c\ , c 2 > 0. 

The PDE for the self-financing hedging strategy takes in this case the following 
form 

1 o 2 S 2 u ss 

2 (1-pciSussY 

This equation was first introduced in ifTOl . The complete description of the 
symmetry group, invariant reductions and invariant solutions to this equation 
as well as to some generalization of this equation are given in (6l, and (HJ. 

2. If the function is given by g(cc) = C2CC C> then the corresponding utility function 
takes the form 

g(a)=c 2 a c \ -> U{x) = -{c 2 xY^ + \. (2.11) 

Under our assumption the utility function acts on (0, +°°) correspondingly the 
function g(ct) = C2Ct Cl should be positive. The assumptions ([T]) - ^ of Defi- 
nition [2j] are satisfied if 

ci,c 2 >0. (2.12) 
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The PDE which corresponds to this function g(tt) has the form 

1 o 2 S 2 u ss u 2 s 

2 [us-ciSu ss Y 

This equation is the special model for hedging strategies in illiquid markets 
which we will study in the section [4] of this paper. 

3. If we take similar to the previous case the function g(cc) in the form g(cc) = 

_ j_ 

C2(p+ka) c < , where a,C2,k G M, k ^ then the corresponding utility func- 
tion is given by 

g(a) = c 2 (p+kay^, -> U{x) = - l -{c 2 x) c ' + (\ + ^y (2.14) 

The utility function acts on (0, +°°) because of that the function g(cc) = C2(p + 
ka)^ 1 ^ 1 should be positive. 



The assumptions of Definition 2. 1 on the the utility function acting on (0, +°°) 
are satisfied if 

c 2 >0, ci€ (-<»,0) U (0,1), k Cl >0. (2.15) 
The PDE for a self-financing trading strategy takes in this case the form 

u t + \ ° 2 ( l+ku s) 2 7 S 2 Uss = o , c, ^ o, (2.16) 
c\ M +ku s +^Su ss 

This equation coincides with some of the models introduced in |fT5l ). This 
model was studied with the Lie group method in [ 4 ] and [ 5 ] . The corresponding 
symmetry algebra admitted by this equation, an optimal system of subalgebras, 
the set of invariant reductions and in some cases also exact invariant solutions 
are presented in these papers. 



3 The Lie algebra admitted by the general model ( 1.7 ) 



To study the symmetry properties of ( |1.7| > we use the method introduced by Sophus 
Lie and developed further in ifTTl . |[T6]l and ITTTl . For the first reading also the books 
|[T9l and can be used which contain many examples. 

The main idea of this method can be formulated in the following way. We in- 
troduce first a jet-bundle of the corresponding order (in this case of the order two, 
because of we have to do with a second order PDE). Then we study all smooth point 
transformations which locally keep the solution subvariety of the studied PDE invari- 
ant. Sophus Lie has proved that instead to look for the symmetry group we can first 
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study the symmetry algebra of the underlying PDE and then with help of an expo- 
nential map we obtain the corresponding symmetry group. We follow these ideas and 
determine the symmetry algebra admitted by the equation. 

Let us introduce now the necessary notations. 
We have two sets of variables, the independent variables t,S which are in the space 
X, (t,S) 6 X and a dependent variable u which belongs to the space U, i.e. u £ U '. 
We consider the second prolongation of the space U. We first introduce the space 
£/(!). It is the Euclidean space endowed with coordinates (u t ,us) which represent all 
first derivatives of u with respect to (t,S). Then we introduce the space Un\, it is the 
space endowed with coordinates (u tt ,u t s,uss) which represent all second derivatives 
of u with respect to (t,S). Now we can define the second prolongation of U and a jet 
bundle of the order two. 

Definition 3.1 The 2—nd prolongation ofU is = U x Um x . 

Definition 3.2 The 2—nd jet bundle of M C X xU, or the jet bundle of order two, 
is =MxUm x Un) = X x U^ 2 \ also called 2—nd prolongation ofM, it is 
also denoted by pA 2 'M. 

We write the PDE < fT77] > in the form 

A(w) =A(x,w (2) ) =A(t,S,u,u t ,u s ,...,u ss ) = 0, wGM (2) , u (1) G U {2) . (3.17) 

Here A is a smooth map from the jet bundle to some Euclidean space R, i.e. 
A : M^ 2 ) — > R. It defines the solution subvariety 

L a = {(i,^):A(j,h( 2 >)=0}cM( 2 ). 

The symmetry group G r of A will be defined by 

G r = {gGDiff(M( 2 ))| g: L A ^L A }, r = dimG r . (3.18) 

The symmetry algebra L r = 3>iff A (M) of the PDE A(w) = can be found as a 
solution of the determining equations 

p A 2 ty(A) = {mod(A(w) = 0)), (3.19) 

where pr^V denotes the second prolongation of an infinitesimal generator V and 
has the following form 

pr^V = Z(S,t,u)j^ + T(S,t,u)-^ + <p(S,t,u)j^ 

+ (p s {S,t,u)^- + (p t (S,t,u)^- (3.20) 
OUs ou t 

+ (p ss (S,t,u)^— + (p s '{S,t,u)^— + (p n (S,t,u)^-. 

du S s ou St du tt 
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Here the last five coefficients can be uniquely defined using the first three (see [16]). 



The determining equations (3. 19 1 is a system of PDEs on the functions £,(S,t,u), 
z(S,t,u) and (p(S,t,u). The solution of this system provides us all infinitesimal gen- 
erators admitted by the studied equation. The set of these infinitesimal generators 
forms a Lie algebra which is called the Lie algebra admitted by the PDE A(w) = 
or the symmetry algebra of this equation. The determining equations and the corre- 
sponding solutions to equation ( |1.7| > are presented in detail in lfT3l . 

We formulate the main result concerning the equation ( |1.7[ ) where we do not put 
any constrains on the form of the function g(tt). We just assume that this function is 
a differentiable one. In this case we obtain the following theorem. 



Theorem 3.1 The equation (1.7 ), where g(cc) is a differentiable function, admits a 



three dimensional Lie algebra L3 spanned by the following generators 

d d d d 

L 3 =<V U V 2 ,V 3 >, V^S^ + u—, V2 = ^, V 3 =3-. (3.21) 

aS au ou at 

The Lie algebra L3 possesses a nonzero commutator relation [V\, V2] = — V2. The Lie 
algebra L3 has a two-dimensional subalgebra L2 =< V\ , V2 > spanned by the gener- 
ators V\ , V2. The algebra L3 is a decomposable Lie algebra and can be represented 
as a semi-direct sum L3 = L2 V3. 



The determining equations ( 3.19| > admit a richer set of solutions if the function 



g(cc) has a special form. We proved (see for details |[T3l ) that equation ( |1.7| admits a 
four dimensional Lie algebra if and only if the function g(ct) has one of the following 
forms: 

g(a) = c 2 e ClU , (3.22) 

g(a) = c 2 a c \ (3.23) 
_ j_ 

g{a) = c 2 (p+ka) c i , ci,c 2 ,& -const. (3.24) 

We discussed in the previous section relations between the given form of the function 
g(cc) and the corresponding utility function. We provided also the PDEs related to 
the special choice of the function g(cc). As we mentioned before the first and the last 
case, i.e. the equations (2. 10 1 and ( |2.16[ ) were studied with the Lie group analysis 



before. The second case, i.e. equation (2.13 1 will be studied in section|4] 



The Lie algebra L3 is the symmetry algebra of equation ( |1.7[ ). Using the usual 
exponential map we can find the symmetry group G3 of this equation as well. It is not 
necessary to find the explicit form of the symmetry group G3 if we will use invariants 
of this group or its subgroups to find reductions and invariant solutions of the studied 
equation. It is enough to know and to use the properties of the symmetry algebra 
which corresponds to the symmetry group. The most interesting reductions we ob- 
tain if we use one-dimensional subalgebras of L3 (and corresponding subgroups of 
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G$). The problem how to choose the optimal system of subalgebras (and correspond- 
ingly subgroups) which give us non-conjugate invariant solutions was solved 1982 
by Ovsiannikov (see [17]). We follow the algorithm of determining of an optimal 
system of subalgebras proposed there. 

We use the notation hf for the subalgebras of the Lie algebra L3 and Hf for the 
corresponding subgroups of the group G3. All three- and four-dimensional solvable 
algebras were classified in [18]. In addition in this paper there are provided optimal 
systems of corresponding subalgebras. We use these results in this and in the next 
section. 



Proposition 3.1 [18]. The optimal system of subalgebras o/Lt, contains the subal- 
gebras shown in Table^ 



Dimension of 
the subalgebra 


Subalgebras 


1 

2 


hf = {V 2 }, h\ = {V lC os(<p) + V 3 sin(<p)}, h 8 3 = {V 2 + eV 3 } 
h g 4 = (V 2 ,V 3 ), h g 5 = (V u V 3 ), h\ = {V l +xV3,V 2 ) 



Table 1: The optimal system of one- and two-dimensional subalgebras of L3, where 
parameters are x € K, < <j> < 71, £ = ±1. 



3.1 The symmetry reductions of the general model ( 1.7 ) 



There are three non-similar subalgebras of the dimension one which we list in Table[2j 
To each of these subalgebras corresponds a one-dimensional subgroup and they form 
the set of non-conjugate subgroups. We use the optimal system of subalgebras for 
invariant reductions of the equation ( |1.7[ ) to ODEs. To obtain reduced form of the 
PDE we take invariant expressions as a new dependent and an independent variables 
and in the new variables we obtain ODEs. In the first line of Table [2] we see that the 
first subgroup Hf describes translations of the dependent variable u and can not be 
used for any reduction. 

In the case of the subgroup Hf the reduction of the given PDE to an ODE is pos- 
sible. But the solutions to the general model ( |1.7[ ) of this type are not very interesting 
because of they have a trivial dependence on time (see the second row, second line in 
Table |2j). The most interesting case we obtain if we use the subgroup Hf. We obtain 



solutions to ( 1.7 ) of the form u(S,t) = SW(lnS — yt) if we are able to solve the first 
order ODE which is represented in the last row and last line of Table [2] It is rather im- 
possible to present a general solution to this equation for an arbitrary function g(a), 
but for some special choices of g(a) it seems to be feasible. 
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Subalgebra 



Invariants 



Transformations/Reductions 




Z = S, W(z)=t 



z = S, W(z) = u-et 



z = InS- yt , 7= cotan(0) 
W{z) = uS- 1 



u — > u + const. 



2e + - 



o 2 z 2 Y' 



(i-z(ing( P Y)yy 

g{pY{z))^cz. 

a 2 (Y'+Y) 



0, Y = W, 



2yY 



0, 



(\-(lng(Y'+Y))') 2 

Y = W, g(Y + Y')^ce z 



Table 2: In the first row we list one-dimensional subalgebras from the optimal system 



of L3 (3.21 1. In the second row we list the corresponding invariants. In the last row 
we provide the transformations or reductions of the PDE ( |1.7[ ) to ODEs. Here e = 
±1, < <j> <7V, a 2 ,p G M + ,c G M are parameters, Y',W' denote the differentiation 
of the corresponding function on the invariant variable z. 



4 Symmetry properties of the special model (2.13) 



We mentioned before that for some special types of the function g(oc) (3.22 1-(3.24 1 



the equation ( |1.7[ ) admits four dimensional Lie algebras. Because two of these cases 
(3.22 >-( 3T24J) were studied before we investigate now the last one, i.e. the symmetry 
properties of the special model (2. 13 1 with the function g(cc) in the form (3.23 ). 
This equation ( 2.13| > contains an arbitrary constant c\ but it does not include any 



more the parameter p. The solution of the determining equation (3. 19 1 defines four 
generators of the Lie algebra L4. We formulate the results in the following Theorem. 



Theorem 4.1 Equation (2.13) admits a four dimensional Lie algebra 

L 4 =<V U V 2 ,V 3 ,V 4 > 



spanned by generators 
Vi = 



s Ts> v > 



v 3 



(4.25) 



± 7 v 4 = l. 

du J du dt 
The Lie algebra L4 possesses the following non-zero commutator relation, \V\ , V3] = 
— V3 . The Lie algebra L4 has a two-dimensional subalgebra L\ =< V\ , V3 > spanned 
by the generators V\ , Vj. The algebra L4 is a decomposable Lie algebra and can be 
represented as a semi-direct sum L3 = L2 V2 V4. 
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To provide non-equivalent reductions of the PDE (2. 13 1 using the Lie algebra 
L 4 we need an optimal system of one-, two- and three-dimensional subalgebras of 
L 4 . Usually just one-dimensional subalgebras give us some interesting non-trivial 
reductions. Also in the case L4 we can use the classification provided in |[T8ll . In the 
paper [18] the Lie algebra L4 ( 4. 1 1 is denoted by L|. 



Proposition 4.1 [18]. The optimal system of the one-dimensional subalgebras con- 
tains four subalgebras, the optimal system of the two-dimensional subalgebras con- 
tains five subalgebras and the optimal system of three-dimensional subalgebras in- 
cludes three subalgebras. All subalgebras from these systems are listed in Table^ 



Dimension of 
the subalgebra 


Subalgebras 


1 


hi = {t/3} , h 2 = {U\ cos0 + £/4sin(0)} , 
h 3 = {U 2 +x(Ui cosO) + t/ 4 sin(0))}, 
h 4 = {U 3 + s(Ui cos(0) + £/ 4 sin(0))} 


2 


h 5 = {U 2 +x(Ui cos(0) + £/ 4 sin(0)),[/ 3 }, 

h(, = {U2 +x(Ui cos(0) + t/4sin(0)), t/i sin(0) — t/ 4 cos(0)} , 

h 7 = {U h U 4 }, 

h = {U3 + e(Ui cos(0) + £/ 4 sin(0)),£/i sin(0) - £/ 4 cos(</>)} , 
h 9 = {Ui,Ui sin(0) — U 4 cos((f>))} 


3 


h w =(U 2 ,U h U 4 ), h 2 u = (U 3 ,U U U 4 ), 

hvi = (U 2 +x(Uicos(^) + U 4 sin((j>)),Ui sin(0) -£/ 4 cos(0),t/ 3 ) 



Table 3: The optimal system of the one-, two- and three-dimensional subalgebras of 
L 4 (4.25 ). Here the parameters are x G R, < <j> < n, £ = ±1. 



4.1 Invariant reductions of the special model (2.131 



Case H\ . The subgroup Hi related to the one-dimensional subalgebra 



does not give rise to any reduction of the PDE (2.13 1 describing the special model of 
the self-financing hedging strategies in illiquid markets. It can be used to modify the 
existing solutions, i.e. if we add an arbitrary constant to a solution it will also be a 
solution of the same equation. 

Case H 2 . The subalgebra h 2 generates after the exponential map a subgroup 
H 2 G G 4 . According to Table[3]it is spanned by the generator 

h 2 = /scos(0)^+sin(</>)^-\, < < %. (4.27) 
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Under the action of the subgroup H 2 expressions 



z = lnS-yt, W{z) 



y = cotan(</>) 6 ' 



(4.28) 



are invariant. 

Remark. For the description of the optimal system it is convenient to use pa- 
rameters in the form cos(0) and sin(</>) to take into account that both coefficients 
are in this way connected to each other. By studying invariant reductions we can 
replace these parameters by 7 = cotan(0) £ 1 as we did in (4.28 1. The both cases 
cos(0) = 0, sin(0) = lead to trivial invariants and reductions and we exclude them 
here and further. 



We use the invariant expressions (4.28) as the new invariant variables. Then we 
obtain from ( |2.13| ) the second order nonlinear ODE of the form 



W' 



(W"-W')W 
(W"-j5W 



a 



i\2 



0, K 



2yc\ 



Cl + 1 



,c h 7€M\{0}. (4.29) 



We excluded both cases c\,y= because these values for the parameters lead to 
trivial invariants and equations. The second term in equation (4.29 1 has a denom- 
inator. We assume that it does not vanish. It means we exclude from the further 
investigations all functions of the type (in the original variables u,S,t) 



u(S,t)=d l S f5 e- 7pt + d 2 , j6 



ci + 1 



7= cotan(</>), d\, d 2 G 



(4.30) 



Now we can multiply both terms of equation (4.29 ) by the denominator and after the 
substitution W = Y ( |4.29| ) reduces to the ODE 

(Y') 2 -Y'Y(2p + K) + Y 2 (p 2 + K) (4.31) 

This is an ODE of the first order and its solution is 

Y(z)=d l e k ^ 1 £ li2 = /3 + |± *Vk(k + 4G3-1)). (4.32) 

The related to ( 4.32| ) solutions of the original equation (2.13 ) take the form 

ii(5,f) =di S kl2 e~ rku2t +d 2 , d h d 2 eR. (4.33) 



We see that the expressions for solutions (4.33) and for the the functions (4.30 1 on 
which the denominator in (2.13 ) does vanish do not coincide identically. 

The solutions (4.33 1 present so called power options or futures. The payoff 
of these options looks similar to usual options, but, for example, instead of the 
usual payoff for a Call option we have at expiry day t = T the following payoff 
Callp ower (S,T) = max[0,S k —B], where B is the exercise price and k is the power 
value. 
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10(1 



Figure 1: Plot of the explicit solutions u(S,t) (4.33) with k\ 
and 1(2 = 1.762 (right panel) for the parameters j8 
(f> = 1.17, c\ = 2.1, d\ = 1, d-2 
and? e (0.1,1). 



1.298 (left panel), 
1.476, k = 0.11, 7 = 0.423, 
0. The variables S,t are in intervals S 6 (0.1, 100) 



Case H3. The third one-dimensional subalgebra in Table [3] is spanned by the 
generator 



h. 



+ x ( Scos(0)— + sin(0) 



dt 



, 0<(j) < 7T,xG M\{0}. (4.34) 



If the coefficient x = 0, then the operator U does not induce any reductions of the 
equation (4.26 1. As the set of invariants of the corresponding group Hi, we take the 
following expressions 

z = InS-yt, W(z) = \nu-St, 8 = (xsin^)) -1 ,/ = cotan(0), y,8 el\{0}. 



We reduce the PDE equation (2.13 ) to a second order ODE using these invariants as 
new variables 



8 - yW' + 



a 2 W' 2 (W" + W a -W) 
2cf (W" + W' 2 -pW) 2 
ci,y,5GK\{0}. 



(4.35) 



We multiply ( 4.35| ) by the denominator of the second term and exclude from the fur- 
ther study all the functions W(z) on which the denominator vanish, i.e. the functions 
(we present this set of functions in the original coordinates u,S,t) 



i{S,t)=d l S^e t{S -^+d 2 e 8t , d u d 2 G 
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After the substitution W' = Y, we reduce equation ( 4.35| > to 

(Y') 2 (8 - yY) + Y'Y (2(5 — yY)(Y — /3) + kY) 
+Y 2 ((8-yY)(Y-P) 2 + kY(Y-1)) =0, 



(4.36) 



2c? 



, q,r,5eM\{0}, 



Solutions of this first order ODE we obtain after the integration 

(5 - yY)dY 



Y [2(8 - yY) (Y - /3) + kY ± ^6Y(Y -Q 
6 = K(4(p-l)y+K), c 



z + di 



(4.37) 



4(j3-l)g K= a^ 
4y(j8 — l)r+ k:' 2cf 



The integrand in (4.37 1 can be transformed into a rational function of a new variable 
T by using the third Euler substitution 



F = 0£(e- T 



2\-l 



(4.38) 



The solution of (4.37 1 contains expressions for roots of a polynomial of a degree 
four and is very space consuming. We do not display this expression here. To ob- 
tain the expression for W(z) and correspondingly for u(S,t), we need one additional 
integration because of the substitution Y = W'. 

Case H4. Let us take the fourth one-dimensional subalgebra listed in Table [3] It 
is spanned by the operator 



/?4 



^ + e r5^cos(0) + ^sin(0) j ) . £ </) £ n.e 11. 



(4.39) 



We find the invariants of the corresponding Lie subgroup H4 if we solve the Lie 
equations, we skip this simple calculations. The invariants in this case H4 can be 
chosen in the form 

z = lnS-yt, W(z) = u-rjt. 



The parameters are 



7 = cotan(0), rj 



sin(0) 



, < < 7T,£ = ±1. 



(4.40) 



Using the invariants (4.40 1 as new invariant variables we obtain the following ODE 
of the second order from (12. 13b 



Tj - yW' + K 



(W"-W')W 2 _ _ 

(W"-pW>) 2 ' K ~2c 2 ' 



(4.41) 
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We exclude from the set of solutions all functions for which the denominator of the 
second term in (437} vanishes, i.e. the family of functions (in old variables u,S,t) 



u(S,t) =c?iS / V 3r? + rjtd 2 , di,d 2 e 



(4.42) 



Then we multiply both terms of equation (4.31 1 with the denominator of the second 
term. After the substitution W' = Y we obtain the first order ODE 

(Y') 2 (ri - yY) + Y'Y(Y(2py+ k) -2]3t]) - Y 2 (Y(p 2 y+ k) - j8 2 T]) = 0. (4.43) 

The set of solutions to equation ( |4.43[ ) coincides with the set of solutions to 

(tj - yY)dY 



VKr((4(j3-i)y+K)r-4(j3-i)rj; 



z + d\ 



(4.44) 



To integrate the expression in ( |4.44| ) we use the substitution (4.38 1, but now we have 
a different expression for £. It means we use the following substitution 



r2\-l 



9 = K(4(P-l)y+K), C 



4(j8-l)T7 
(4(J8-I)y+K) 



(4.45) 



After the substitution (4.45 1 equation (4.44 1 takes the form 



rdr 



-4t] 

~dTJ b 2 T 2 ±T + b 



+ 4y 



TdT 



{d-T 2 )(b 2 T 2 ±T + b Q ) 

-1 



z + di 



b 2 = 2p(6Q- L , b = K + 2py-2P£- 1 , d x € 

and can be easily integrated. 

We display the solution to ( |4.44| ) with the parameter e = 1 



(4.46) 



(4.47) 



2p(p 2 y+K)z + d 1 =2(p 2 y+K)lnY 
-J8 Vd In - + v / 0F(0F-ai)^ 

+(J8 -2)Kln (F( j 8 2 ( J c + 2 7 ( j 8 - 1))) - 2j3 2 (j3 - 1)t] +j3(j3 -2) V / P^) , 



where ai = 4fc(/3 - 1)tj, = 4jc4(j3 - l)y+ jc 2 ), G E. The solution to ( |4.44 
with £ = — 1 has a similar form and we skip this expression. The plot of a solution to 
(4.44 1 for e = 1 is given in Figure [2] 

To obtain a solution to equation ( |4.41 1 we should first invert the function z = 
f(Y) of the type (4.47) and then integrate W(z) = jY(z)dz. In the general case 
neither the expression (4.47 ) nor the similar expression with e = — 1 can be inverted 
explicitly. But we have possibility to solve this problem numerically and in this way 
solve equation (|2.13 1. 
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5 Conclusion 



We studied with the Lie group analysis the general equation ( |1.7[ ) and obtain the 
three dimensional symmetry algebra L3 (|3.21|). Using the Lie algebra L3 and the 



corresponding optimal system of subalgebras we provide the complete set of non- 
equivalent reductions of the PDE ( |1.7[ ) to ODEs. The reductions are listed in Table[2j 
In addition we were able to figure out all possible forms of the function g(cc) 
for which equation ( |1.7[ ) admits an extension of the symmetry algebra L3. All such 
functions are listed in ( |3.22[ >-( 3.24 1. From the three equations with the special form of 
the function g(cc) two equations were deduced earlier by different authors iTTOl . |[T5l 
and studied with the Lie group analysis in EJ-lSl. The new equation in this list is 



equation (2.13 1 which corresponds to the function g(a) (3.23 1. The main properties 
of the symmetry algebra L4 for this equation are represented in Theorem |4.1| The 



invariant reductions of (2.13) are described and studied in section [4] Using, similar 
to the general case, an optimal system of subalgebras we describe the complete set 
of non-equivalent reductions of the PDE to different ODEs. In all these cases we are 
able to solve the ODEs and obtain solutions in exact form. We skipped in one case 
the exact formula because it is too voluminous. 



It is interesting to remark that the models (2.10l, (2.16l and (2.13) where intro- 



duced under different finance-economical assumptions. Now we clarified how deep 
they are connected from mathematical point of view. Other interesting remark con- 
cerns the structure of the corresponding Lie algebras. The symmetry Lie algebras of 



the models (2.10l, (2.16l and (2.13) are isomorphic. It is easy to understand because 



of all of them are solutions to the same determining system (3.19). Less evident is 



why the Lie algebra admitted by the risk-adjusted pricing methodology model stud- 
ied in [3] is also isomorphic to the previous ones. It seems that many of the pricing 
options models in illiquid markets posses similar algebraic properties. 
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